Astronomy & Astrophysics manuscript no. AA200912599 


© ESO 2009 


November 3, 2009 





OS 
o 
o 

(N 
> 

o 



X 

Oh 
6 

(N 
> 
00 
On 



INTEGRAL probes the morphology of the Crab nebula in hard 

X-rays/soft 7-rays 

D. Eckert 1 ' 2 , V. Savchenko 1 , N. Produit 1 & C. Ferrigno 1 - 3 

'ISDC Data Centre for Astrophysics, Geneva Observatory, University of Geneva, 16, ch. d'Ecogia, CH-1290 Versoix, Switzerland 
2 INAF/IASF Milano, Via E. Bassini 15, 20133 Milano, Italy 

3 IAAT, Abt. Astronomie, Universitat Tubingen, Sand 1, 72076 Tubingen, Germany 



Received May 29, 2009/ Accepted October 08, 2009 



ABSTRACT 



Aims. We use the IBIS/ISGRI telescope on-board INTEGRAL to measure the position of the centroid of the 20-200 keV emission 
from the Crab region. 

Methods. We find that the astrometry of the IBIS telescope is affected by the temperature of the IBIS mask during the observation. 
After correcting for this effect, we show that the systematic errors on the astrometry of the telescope are of the order of 0.5 arcsec. 
In the case of the Crab nebula and several other bright sources, the very large number of photons renders the level of statistical 
uncertainty in the centroid smaller or comparable to this value. 

Results. We find that the centroid of the Crab nebula in hard X-rays (20-40 keV) is shifted by 8.0 arcsec with respect to the Crab 
pulsar in the direction of the X-ray centroid of the nebula. A similar shift is also found at higher energies (40-100 and 100-200 keV). 
We observe a trend of decreasing shift with energy, which can be explained by an increase in the pulsed fraction. To differentiate 
between the contribution of the pulsar and the nebula, we divide our data into an on-pulse and off-pulse sample. Surprisingly, the 
nebular emission (i.e., off-pulse) is located significantly away from the X-ray centroid of the nebula. 

Conclusions. In all 3 energy bands (20-40, 40-100 and 100-200 keV), we find that the centroid of the nebula is significantly offset 
from the predicted position. We interpret this shift in terms of a cut-off in the electron spectrum in the outer regions of the nebula, 
which is probably the origin of the observed spectral break around 100 keV. From a simple spherically-symmetric model for the 
nebula, we estimate that the electrons in the external regions of the torus (d ~ 0.35 pc from the pulsar) reach a maximal energy 
slightly below 10 14 eV. 

Key words, pulsars: individual: Crab - ISM: supernova remnants - Gamma rays: observations 
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The Crab nebula is the brightest astrophysical source in the y- 
ray sky (E > 30 keV). It is the pro t otypic al pulsar-wind nebula 
(PWN, see e.g., iGaensler & Sland d2006l) for a review), where 
material ejected from the central pulsar at different epochs in- 
teracts, producing st rong shocks that accele rate electrons up to 
energies ~ 10 15 eV dKennel & Coronitilll984j) . The synchrotron 
emission produced by the large populat ion of non-thermal el ec- 
trons in a magnetic field B ~ 0.3 mG dMarsden et al.llT984b is 
detected over more than 10 orders of mag nitude in the electro- 
magnetic spectrum from radio to y- rays dAtovan & Aharonianl 
1996: IVolpi et al.ll2008l: [Hester1l2008h . 

In soft X-rays, th e size of the nebula is ~2 arcmin 
dWeisskopf et alj|2000l) . It consists of a characteristic jet+torus 
structure, probably corresponding to relativistic outflows along 
the rotation axis and the equator of the pulsar. High-resolution 
Chandra observations also revealed the presence of an inner 
ring, probably associated with the conversion of the relativistic 
pulsar wind into a synchrotron-emitting plasma. The measured 
photon index decreases with radius from T ~ 1.6 around the 
pulsar down to T ~ 3.0 in the outer p arts of the X-ray image 
dMori et al.ll2004IWillingale et alj|2001l) . which implies that the 
most energetic electrons are constantly injected by the pulsar in 
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the inner regions of the nebula. The steeper photon index in the 
outer regions can be explained by synchrotron cooling. This in- 
terpretation is supported by the observation of rapid X-ray vari- 
ability in the inner ring dHester et al J 12002). which is probably 
due to the presence of strong quasi-stationary shocks. The pul- 
sar itself exhibits a hard spectrum T ~ 1 .6 and a period P ~ 30 
msec with a double-peak profile. 

In the hard X-ray/soft y-ray range, early HEAP A-4 observa- 
tions inferred a photon index of F — 2.15 (Jung 1989). Around 
100 keV a spectral break was detected, and above this energy 
a photon index of ~2.5 was found. This measurement agrees 
with the stee per spectral indice s measured at higher energies b y 
COMPTEL dStrone et alJ 19931) and EGRET dNolan et alJl993l) . 
Therefore, observations of the Crab complex around the break 
energy by modern experiments are important to constrain parti- 
cle acceleration models. Recently, a polarized y-ray s ignal from 
the Crab was mea sured by the SPI dDean et alj |2008) and IBIS 
dForot et alJl2008l) instruments on INTEGRAL. The polarization 
was found to be co-aligned with the spin orientation, thus in- 
dicating a possible association of the polarized signal with the 
inner jet. The broad-band coverage of INTEGRAL also permit- 
ted a st udy of the evolutio n of the pulse profiles as a function of 
energy dMineo et al.l2006l) . A phase-revolved analysis detected a 
significantly harder signal (F ~ 1.8) during the interpulse phase 
compared to the peaks (F ~ 2.2). However, the results presented 
there only considered the emission from the pulsar, and no in- 
formation was given about the spectrum of the nebula. 
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In this paper, w e use the coded-apert ure IBIS telescope 
dUbertini et al.l2003l) on-board INTEGRAL dWinkler et al.l2 003) 
to measure the position of the hard X-ray/soft y-ray centroid of 
the Crab pulsar/nebula complex. In Sect. [2] we analyse the as- 
trometry of the IBIS telescope and show that the point-source 
localization accuracy of the instrument depends on the temper- 
ature of the mask. We correct for this effect and assess the level 
of systematic uncertainties in the astrometry of IBIS. In Sect. 
[3] we report precise measurements of the position of the Crab 
pulsar/nebula and present our results. The implications of these 
findings are discussed in Sect. [4] 

2. Astrometry of the IBIS telescope 

2. 1 . IBIS misalignment correction 

The IBIS te lescope, and in p articular its low-energy detector 
layer ISGRI dLebrun et al.ll2.Q03D , operates in the hard X-ray/soft 
y-ray range ( 15-400 keV) and has an angular resolution of 12 ar- 
cmin FWHM. While this is much larger than the angular extent 
of the Crab nebula, the point source location accuracy of the in- 
strument scales with the inverse of the detection significance of 
the source 

FWHM 



where R = 2.43 is the r atio of the size of mask to detector pixels 
dGoldwurmetalJl200ll) . so in principle for a very bright source 
such as the Crab nebula (S/N ~ 3, 000 in the 20-40 keV band) 
IBIS/ISGRI can measure the position of the peak with an accu- 
racy of ~ 0.1 arcsec . However, it was noticed early in the mission 
(IWalteret al.l l2003i that the axis of the star trackers was offset 
with respect to the spacecraft axis, introducing a systematic er- 
ror in the positions measured by the telescope. To correct for this 
effect, a misalignment matrix performing a rotation between the 
star tracker and spacecraft axes was introduced, which allowed 
for the localization of astrophysical sources with a reasonable 
accuracy. However, for very bright sources (e.g., Crab, Cyg X-l, 
Sco X-l) a systematic shift of ~ 7 arcsec with respect to the true 
source position was present, and thus for these sources it was 
not possible to measure the position with an accuracy greater 
than this value. 

In the conte xt of the development o f the Offline Scientific 
Analysis (OSA, ICourvoisier et al.ll2.Q03l) version 8.0, we inves- 
tigated this effect analyzing a large sample of Cyg X-l obser- 
vations (~ 900 individual pointings (Science Windows, ScWs) 
each of typical duration 2-3 ksec) and computed the alignment 
of the instrument from scratch. In particular, we investigated the 
effects of the temperature of the IBIS mask on the source lo- 
calization accuracy. Indeed, the dilatation of the IBIS mask and 
its supporting structure produced by the relative position of the 
Sun might have an influence on the alignment of the instrument. 
Figure Q] shows the relative declination of Cyg X-l with respect 
to the catalog position (in arcsec) as a function of the mask tem- 
perature. A linear trend can clearly be seen in the figure, and 
hence it seems that the mask temperature has a significant effect 
on the source localization accuracy of IBIS/ISGRI. 

To correct for this effect, we introduced a linear dependence 
of the three rotation angles (i.e., Euler angles around the three 
spacecraft axes) on the mask temperature, and fitted the an- 
gles using a least squares method to achieve the closest possi- 
ble agreement between the observed and catalog positions. We 
found that only one angle depends significantly on temperature, 
and therefore only one additional parameter is needed to correct 




Mask temperature [C] 



Fig. 1. Relative Declination of Cyg X-l with respect to the cata- 
log position (in arcsec) as a function of the IBIS mask tempera- 
ture (in °C). The black line shows a linear fit to the data. 



for the temperature effect. The dependence of this angle on the 
temperature was measured accurately by the fit. Finally, we im- 
plemented these results into the scientific analysis software. 

2.2. Test of the IBIS astrometry on several point sources 

To validate our method and test it on well-known cases, we se- 
lected a large sample of data comprising 2,017 individual point- 
ings and containing at least one of the 3 brightest point sources in 
the ISGRI band (Sco X-l, GRS 1915+105, and V0332+53 dur- 
ing its 2005 outburst, but excluding Cyg X-l since it was used 
as calibration source) and measured the position of the sources 
in each ScW (20-40 keV band). 

We restricted our observations to pointings with low off-axis 
angle (< 8°), since the width of the IB IS PSF increases signifi- 
cantly for larger angles dGros et al. I I2003I Using these data, we 
constructed histograms of the relative shift in RA and Dec of the 
measured source position with respect to the catalog one. The re- 
sulting plots are presented in Fig. [2] Both distributions are well- 
represented by a Gaussian profile with a similar standard devia- 
tion cr = 8.2 + 0.1 arcsec (RA) and <x = 8.1 + 0.1 arcsec (Dec). 
The mean of the distribution is found to be fi = -0.35 + 0.21 
arcsec (RA) and // = 0.51 + 0.23 arcsec (Dec). Similar shifts (al- 
though less accurate statistically) are obtained when each source 
is analysed individually. Based on these results, from now on we 
adopt the value of 0.5 arcsec as an estimate of the systematic 
uncertainties in the astrometry of the IBIS telescope. 

The reader might argue that the systematic errors adopted 
here are smaller than the uncert ainties in the attitu de reconstruc- 
tion of the telescope (3 arcsec. IWalter et al]|2003l) . However, it 
should be noted that the uncertainties in the attitude reconstruc- 
tion affect each pointing individually in a non-preferential way. 
Given that we are analysing a large sample of pointings, this in- 
troduces a hard limit to the width of the distribution (~ 8" » 3"), 
but does not affect the mean of the distribution. 

In addition, one might also wonder whether the alignment 
of the instrument depends on energy, since the alignment was 
performed using low-energy data (20-40 keV) where count rate 
and signal-to-noise ratio are maximal. However, we note that 
the misalignment between the instrument and the star trackers 
is caused only by geometrical and thermal effects, which do not 
depend on the energy of the incoming photons. Therefore, for 
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Fig. 2. Test of the IBIS/ISGRI astrometry on 3 of the brightest 
point sources in the 20-40 keV band (Sco X-l, GRS 1915+105 
and V0332+53). The plots show the relative shift in RA (top) 
and dec (bottom) measured by the telescope with respect to the 
catalog position. Both histograms are fitted by a Gaussian profile 
(black curve). 



the remainder of the paper we use the same alignment and sys- 
tematic errors in all energy bands. 



3. Results 

3.1. Measuring the centroid of the Crab nebula 



As shown in Sect. l2.2l we are now confident that we can measure 
the position of sources with a level of systematic uncertainties 
of about 0.5 arcsec. Since the pulsed fraction is less than 20% 
in the INTEGRAL energy range, while it is ~10% in the 1-10 
keV, the hard X-ray emission is dominated by the nebula and 
thus IBIS/ISGRI can measure its position with excellent accu- 
racy. We selected a data set from the ISDC archival] only for 
pointings for which the offset angle of the Crab with respect 
to the spacecraft axis was at most 8.0 degrees, and we filtered 
out the observations affected by Solar flares. Overall, our sam- 
ple comprises 1,141 individual ScWs of the Crab nebula, of a 
total effective exposure time of ~ 1.8 Msec. 

The analysis was performed with the standard OSA software 
v.8.0, including the correction of thermal effects presented in 
Sect. 12.11 For each pointing, an image was created in the 20-40, 




o 5 10 
Relative R.A. [arcsec] 



http : //isdc . unige . ch/?Data+browse 



Fig. 3. Relative right ascension (R.A.) of the best-fit position of 
Crab in the 20-40 keV band with respect to the position of the 
Crab pulsar. The black solid line shows a Gaussian fit to the data. 



40-100, and 100-200 keV bands, and the position of the Crab 
nebula was fitted by a 2-dimensional Gaussian with the width 
of the IBIS PSF. The best-fit position was then collected, and 
histograms were computed. Figure [3] shows the distribution of 
the Crab best-fit position (in RA) in the 20-40 keV band with 
respect to the Crab pulsar. The distribution is obviously shifted 
from the Crab pulsar by ~5.5 arcsec to the west. A similar shift 
is also found in declination by ADec~6.3 arcsec towards the 
north. We then fitted the distributions by a Gaussian profile to 
extract the best-fit positions (see Table Q]). In the 20-40 keV 
band, the centroid of the Crab nebula is found to be at RA = 
83.6317 ± O.l'W + 0.5" syst , dec = 22.0162 + 0.1" stat + 0.5" syst . 
This corresponds to an angular distance to the Crab pulsar of 
8.0 + 0.1 s tat ± 0.5 sys t arcsec. The fitted position is away from the 
Crab pulsar by more than 13o\ 

Given that t he timing resolution of ISGRI (90yus, 
iLebrun et al.ll2003l) is sufficient to reso lve the on- and off- pulse 
emission from the pulsar (P = 30 ms, iTennant et alJl200lT) . we 
extracted phase-resolved images of the Crab nebula for both 
the on- and off-pulse phases (see Table Q]). For convenience, 
the 2 phase bins were chosen to be equally-long, th e on-pulse 
bin co mprising both peaks and the interpulse (see Mine o et"aTI 
2006). More specifically, the on-pulse bin includes phases 0.95- 
0.45 (see Fig.[5j. In the 40-100 and 100-200 keV bands, we were 
also able to measure the centroid of the total emission with suf- 
ficient accuracy (see TableQ]). 



3.2. Comparison with soft X-rays 

To visualize the best-fit positions presented in Sect. 13. 11 we used 
an archival high-resolution Chandra observation, and extracted 
images in the full band (0.5-10.0 keV) as well as in a soft (0.5- 
2.0 keV) and a hard band (2.0-10.0 keV). In the Chandra obser- 
vation, the very strong pile-up effect at the position of the pul - 
sar leads to the rejection of the events ( Weissko pf et al.l l2000). 
which allows us to analyse the nebula individually. To estimate 
the centroid of the nebular emission, we convolved the full-band 
Chandra image with a very broad Gaussian and fitted the result- 
ing image with a 2-dimensional Gaussian. The resulting best- 
fit position for the centroid lies 17.3" away from the pulsar, in 
the torus and along the axis of the jet. Figure [4] shows the best- 
fit positions with 1 <x error circles shown in Tab. Q] for the total 
emission and the on- and off-pulse cases, superimposed on the 
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Table 1. IBIS/ISGRI best-fit positions of the Crab in 3 energy 
bands. The phase-resolved positions (on- and off- pulse) are also 
reported. The lcr error on the position (3 column) is the sum of 
the statistical and systematic uncertainties, where the systematic 
error is fixed to be 0.5" (see Sect. 12. 2\ . The last column gives the 
angular distance between the measured position and the pulsar 
(in arcsec). 



Chandra image of the region. We can see that all the measured 
positions lie less than lcr away from the line joining the Crab 
pulsar to the X-ray centroid of the nebula, which corresponds 
roughly to the jet line from the pulsar. This indicates that the 
measured positions are a combination of the pulsar position with 
the position of the centroid of the nebula. 

Interestingly, we can see in Fig.[4]that all the measurements 
are found to be significantly offset from the centroid of the X- 
ray emission. Given that the pulsed emission accounts only for 
~ 15% of the flux, this indicates that the centroid of the emission 
from the nebula, i.e., its morphology, varies relative to the soft X- 
ray band. This statement is confirmed by the measurement of the 
off-pulse emission (where the signal is completely dominated by 
the nebula), which is also found to be significantly offset from 
the X-ray centroid of the nebula. 

3.3. Evolution of the pulsed fraction 

To compare the results described above with the expected val- 
ues, we investigated the evolution of the pulsed fraction with 
energy. Indeed, given that the measured centroid is the sum of 
the point-like signal from the pulsar and of the unresolved neb- 
ular emission, any evolution of the pulsed fraction with energy 
will influence the position of the centroid. 

We perfo rmed a phase-r esolved analysis using the 
ep hemeris of | Lyne et al.l (fl9 93lPl and exploiting the method 
of [Segreto & Ferrigno (20071). We collected all the Crab ob- 
servations within an off-axis angle of 8 degrees (exposure of 
about 1.8 Ms) and produced a folded light curve with 200 phase 
bins, using energy bins of width 0.4787 keV. Figure [5] shows 
the resulting background-subtracted pulse profiles in the 20-40, 
40-100, 100-200, and 200-500 keV bands. The characteristic 
double-peak profile is clearly seen. In addition, the figure shows 
that during the off-pulse phase (0.45-0.95) the measured flux is 
nearly constant, which indicates that the pulsar contributes very 
little to the emission during the off-pulse phase. 

We then determined the pulsed fraction as a function of en- 
ergy following the definition (F on - F ff)/(F on + F q), where 
F on is the total fluence during the pulse phase bins (phase 0.95- 
0.45), and F q is the fluence during the off-pulse bins (phase 0.6- 
0.9). In the energy range 20-200 keV, we found that the pulsed 
fraction depends linearly on log E with slope 0.055 + 0.001 (see 
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Fig. 4. Chandra/ ACIS-S image of the central region of the Crab 
nebula, with the best-fit ISGRI positions in the 20-40 (band 1), 
40-100 (2) and 100-200 keV (3) bands overlayed. The red cir- 
cles show the best-fit positions and lcr error circles of the total 
emission, while the black (yellow) circles show the measured 
positions of the off-pulse (on-pulse) phase bins. The position of 
the Crab pulsar and of the X-ray centroid of the nebula are also 
displayed. The size of the blue arrow is 17.3". For comparison, 
the inset shows a global view of the nebula. The black rectangle 
represents the size of the region displayed here. 



Fig Et- The pulsed fraction increases significantly from 0.13 at 
20 keV up to 0.18 at 150 keV, i.e., by about 30%. 

To relate the measured pulsed fraction to the predicted posi- 
tion of the centroid, we simulated 10 6 photons emanating from a 
point source at the position of the pulsar and an extended source 
with the extension of the Crab nebula, with a centroid shifted by 
17.3" as found from Chandra data. We convolved the resulting 
photon distribution with the large beam of ISGRI and fitted the 
distribution with a Gaussian to determine the position of the cen- 
troid. Unsurprisingly, our simulations indicate that the expected 
centroid depends linearly on the pulsed fraction. Therefore, us- 
ing Fig. [6] the pulsed fraction can be directly compared to the 
observed shift. 



3.4. Relative positions 

In addition to the measurement of the centroid independently in 
the different bands and pulse phases, it is also possible to mea- 
sure the difference between several parameters in each point- 
ing individually. In this case, the systematic error of 0.5 arcsec, 
which probably originates in uncertainties in the spacecraft at- 
titude reconstruction, cancels out, so it is more accurate to cal- 
culate the difference between 2 quantities in each pointing and 
calculate the mean rather than compute the mean of each 2 quan- 
tities and then their difference. This allows us to measure the 
difference between the positions during the on- and off-pulse 
phases with very good accuracy. Moreover, we have seen in Sect. 
13.21 that all the positions lie along the jet line, so it is conve- 
nient to compute the difference between 2 quantities along this 
line. Figure [7] shows the distribution of the difference between 
the on- and off-pulse emission in the 20-40 keV band for each 
ScW, measured along the jet line. Fitting the distribution with a 
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Fig. 5. Background-subtracted pulse profiles of the Crab in the 
20-40, 40-100, 100-200 and 200-500 keV bands. In the bottom 
left corner, the dashed line shows the off-pulse phase bin (0.45- 
0.95) while the solid line shows the on-pulse phase bin. 
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Fig. 7. Distribution of the difference between the on- and off- 
pulse position measurements in each individual pointing, mea- 
sured along the jet line. The distribution is significantly shifted 
by 2.28 + 0.15 arcsec. 



Quantities measured 



Relative angular distance ["] 



On-pulse/Off-pulse, 20-40 keV 
On-pulse/Off-pulse, 40-100 keV 
On-pulse/Off-pulse, 100-200 keV 
Total, 40-100/20-40 keV 
Total, 100-200/20-40 keV 
On-pulse, 40-100/20-40 keV 
On-pulse, 100-200/20-40 keV 
Off-pulse, 40-100/20-40 keV 
Off-pulse, 100-200/20-40 keV 



2.28 : 

2.58: 
3.65: 
1.50: 
2.14: 
1.25 : 
2.84: 
1.05 : 
1.51 : 



0.15 
0.19 
0.98 
0.15 
0.49 
0.19 
0.63 
0.24 
0.84 



Energy [keV] 



Table 2. Mean of the relative distance between two quantities, 
measured along the jet line (see text). 



We can also see from Table [2] that both the on- and off-pulse 
emission shift towards the pulsar with energy. While this is per- 
fectly expected for the on-pulse emission as the signal from the 
pulsar becomes stronger, this argument cannot be valid for the 
off-pulse emission, where by construction the pulsed emission is 
suppressed. Therefore, the relatively significant shift (4.3<x) ob- 
served between the off-pulse emission in the 20-40 and 40-100 
keV bands must be caused by a change in the morphology of the 
nebula. 



Fig. 6. Pulsed fraction of the Crab pulsar/nebula complex as a 
function of energy. The solid line shows a linear fit to the data, 
with lcr error given by the dashed lines. 



Gaussian, we find a very significant shift (15<x) of 2.28 + 0.15 
arcsec between the on- and off-pulse phases. The standard devi- 
ation of the distribution is found to be smaller than that of the 
2 distributions individually, which proves that the systematic ef- 
fects were at least partially cancelled. Perpendicular to the jet 
line, no shift is found (A = 0.16 ± 0.15). This confirms that 
the centroid is moving along the jet line. A similar analysis can 
also be performed for the positions measured in different energy 
bands. All the measured relative shifts are presented in Table[2] 
From Table |2l we can see that the difference between the 
on- and off-pulse phases seems to increase with energy. Since 
the pulsed fraction increases by ~ 30% (see Fig. |6]), this ef- 
fect can be easily explained by the increasing pulsed emission. 



4. Discussion 

As it has been shown above, IBIS/ISGRI is capable of measuring 
the position of the X-ray centroid of the Crab pulsar/nebula with 
unprecedented accuracy (< 1") in the 20-200 keV band, for the 
total emission as well as for the phase-resolved cases. The best- 
fit positions all lie along a line joining the X-ray centroid of the 
nebula to the Crab pulsar. 

To obtain physical insight using these results, we used 
Chandra images in the soft (0.5-2 keV) and hard (2-10 keV) 
bands to create a hardness ratio map (2-10/0.5-2). We ex- 
tracted Chandra/ACIS RMFs and ARFs for this observation 
and used these spectral responses to convert between the 2- 
10/0.5-2 hardness ratio and the corresponding photon index. We 
fixed the galactic absorption in the direction of the source to 
the value measured by XMM-Newton (n H = 3.45 x 10 21 cirT 2 , 
IWillingale etafl d200ll) . in agreement with the 21cm value) and 
simulated Chandra spectra for a wide range of photon indices 
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(1.4-3.0) to compute the corresponding hardness ratio. We then 
fitted the simulated data points using a 3 rd order polynomial 
and used the best-fit function to convert the hardness ratio map 
into a phot on-index map. Our resulting photon-index map agrees 
with that of M ori et al.l ^2004). Finally, assuming pure power-law 
spectra we used our photon index map together with the hard 
band image to extrapolate to higher energies. As a result, we ob- 
tained extrapolated high-resolution maps in the 20-40, 40-100, 
and 100-200 keV bands, for which we measured the centroid of 
the nebula. The centroid of the extrapolated maps can be com- 
pared directly to the off-pulse emission measured by IBIS/ISGRI 
and presented in Tab. Q] To estimate the uncertainty in the ex- 
trapolated centroids, we simulated 1,000 extrapolated maps with 
soft X-ray flux and photon index normally distributed in each 
pixel, and we computed the distribution of the simulated cen- 
troids. The lcr error is then given by the standard deviation of 
the distribution. The resulting lcr error is very small, from 0.06" 
in the 20-40 keV band up to 0. 15" in the 100-200 keV band. 

Comparing the extrapolated centroids with the ISGRI off- 
pulse measurements, we find that the emission is significantly 
offset. In the 20-40 keV band, the ISGRI off-pulse emission is 
measured to be 4.5" away from the extrapolated position (6 Act), 
in the direction of the pulsar. A similar difference (4.8", 6<x) is 
found in the 40-100 keV band. Consistently, we have seen that 
the centroid of the off-pulse emission shifts towards the pulsar 
with increasing energy (see Sect. 13.4b . This indicates that ISGRI 
probes different physical regimes than Chandra. Only 2 different 
explanations of this effect can be invoked: either the centroid of 
the nebula is moving towards the pulsar, or the steady flux from 
the pulsar increases dramatically. 

According to lTennant et al.l d200ll) . the minimum soft X-ray 
flux of the pulsar is < 1% of the peak flux. Moreover, the authors 
conclude that the origin of the emission is probably non-thermal, 
and is due to the same emission mechanism as the pulse, so the 
steady flux should vary in a similar way as the peak flux. To pro- 
duce a shift of 4.5" in the 20-40 keV band compared to the pre- 
dicted position, the steady flux from the pulsar should represent 
~ 35% of the total flux, i.e. around 3 orders of magnitude greater 
than the value measured by Chandra in soft X-rays. Therefore, 
we can reasonably exclude the possibility that the detected shift 
would be caused by an increase in the persistent flux from the 
pulsar. As a result, we conclude that the shift between the pre- 
dicted and the measured position of the Crab nebula is produced 
by a change in morphology of the nebula. 

In the Crab nebula, the most energetic electrons are ex- 
pected to be produced at the w ind termination shock, in the in- 
ner ring ( Weisskopf et al. 2000). With increasing distance from 
the termination shock, the cut-off energy in the electron spec- 
trum decreases because of synchrotron cooling (see Fig. 3 of 
lAtovan & Aharonianll 19961) . In X-rays, the nebula is ~ 2 times 
smaller than in the optical, and therefore the synchrotron emis- 
sion in the outer regions probably cuts off in the optical/UV. 
Following this picture, it is natural to expect that a similar effect 
is happening when we compare the hard X-ray/soft y-ray band 
with the soft X-rays. In the 20-40 keV band, we see in Fig.|4]that 
the centroid of the emission is found to be very close to the inner 
ring, which corresponds to the region where the most energetic 
electrons are produced. This indicates that in the outer regions 
of the X-ray image the spectrum cuts off in the X-ray/hard X- 
ray range. This interpretation is supported by the presence of the 
well-known sp ectral bre ak at ~ 100 keV in the total spectrum 
of the source (Juna fl989l) . which corresponds to an energy of 
~ 10 14 eV for the electrons. In this framework, the measure- 
ments of the centroid of the emission around the break energy 
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Fig. 8. Relative shift of the centroid of the nebula with respect 
to the pulsar position in the 20-40 keV band as a function of the 
cut-off radius R cu , (see text). The black lines show the acceptable 
range (lcr) for our measurement of the off-pulse centroid in the 
20-40 keV band. 

with arcsecond accuracy imply that the spectral break is due to 
a cut-off in the electron distribution away from the termination 
shock. 

To validate this hypothesis, we extrapolated the Chandra im- 
age to higher energies using this model rather than a pure power 
law. For simplicity, we assumed that the source is spherically 
symmetric around the pulsar and defined a radius R cut from the 
pulsar such that for R < R cu , the emission is described by a pure 
power law, while for R > R cul an exponential cut-off is intro- 
duced. For the cut-off energy, we assumed a linearly-decreasing 
profile from E cut = 40 keV at R = R cut (the maximum of the 
corresponding ISGRI band) down to E cu , = 8 keV in the out- 
skirts (the upper limit of the Chandra response). In physical 
terms, the characteristic synchrotron frequency of electrons with 
a Lo rentz factor y in a magnetic field B ~ 0.3 mG is given by 
(e.g. lBlumenthal & Gouldl (fl97Qh ^ 

V -~ 7X1 ^0tLg) (2) 

so we study electrons with a maximal energy E ~ 10 13 — 10 14 eV. 
We then used this model to extrapolate the 2-10 keV Chandra 
image to the 20-40 keV band and computed the centroid of the 
extrapolated image for a range of values of R cut . Figure [8] shows 
the predicted relative position of the centroid of the nebula with 
respect to the pulsar position, as a function of the cut-off ra- 
dius R cut . Surprisingly, for very small values of R cu , the shift de- 
creases with radius. This can be explained by the fact that for 
very small R cut a cut-off is present throughout the entire nebula, 
and therefore the morphology is not changing much. The func- 
tion reaches a minimum at R cul ~ 12 arcsec, which corresponds 
to the radius at which the contrast between the central region 
(power law) and the outer parts (cut-off power law) is the great- 
est. For larger radius, the shift again increases as expected. 

In the 20-40 keV band, the off-pulse emission is shifted by 
~9 arcsec compared to the pulsar (see Tab. [TJ. From Fig. [8] we 
see that 2 different values of R cu , match with value, i.e., R cut ~ 5 
arcsec and R cu , ~ 25 arcsec. However, the smaller value corre- 
sponds to the very inner region of the nebula, and therefore a cut- 
off at this radius would have strong implications on the spectral 
shape, which is not seen, so we can safely reject the first solu- 
tion. As a result, we find that the cut-off radius that corresponds 
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to the observed centroid is R cut ~ 25 arcsec . At the distance of 
the Crab nebula (d ~ 2 kpc, iTrimbld d!968l) ). this corresponds 
to a radius of ~ 0.35 pc, i.e., to the outer regions of the torus. 
This proves that a cut-off effect in the outer regions of the neb- 
ula can reproduce the observed shift. Our modelling implies that 
at a radius R ~ R alt , the synchrotron cut-off energy is E ~ 40 
keV, which corresponds to a maximal energy for the electrons of 
~ 6 x 10 13 eV. 

Of course, these results are strongly model-dependent and 
our model is overly simplistic, in particular because of the as- 
sumption of spherical symmetry. Moreover, the output of our 
model depends on several input parameters, in particular on the 
choice of the maximum and minimum cut-off energy, so the de- 
rived values should be considered with care. However, the result- 
ing numbers point towards a reasonable scenario for the electron 
distribution in the nebula. Indeed, our results are in excellent 
agreement with detail ed theoretical modelling of the source us - 
ing MHD simulations dVolpi et alJl2008tlDel Zanna et alj|2006l) . 
which predict that the jet and the external regions of the torus 
should disappear in hard X-rays/soft y-rays (see Fig. 4 of 
Volpi et all l2008f) . Therefore, our measurements qualitatively 
confirm the predictions of MHD simulations. 

5. Conclusion 

We have used the IBIS/ISGRI instrument on-board INTEGRAL 
to measure the centroid of the Crab pulsar/nebula complex in 
hard X-rays/soft y-rays with unprecedented accuracy, with the 
aim of investigating the behaviour of the relativistic electron 
population around the break energy (.Ebreak ~ 100 keV). Based 
on our understanding of the dependence of the astrometry of 
IBIS on the temperature of the mask, we have shown in Sect. 12.21 
that despite its moderate angular resolution (12 arcmin FWHM), 
in the case of sources detected with high S /N IBIS/ISGRI can 
measure the position of astrophysical sources with an accuracy 
~ 0.5 arcsec. Applying this method to the Crab, we found that 
the emission is offset significantly with respect to the position 
of the pulsar, by ~ 8" in the 20-40 keV band. We were also 
able to measure the position of the on- and off-pulse phases in- 
dependently. All the measured positions are found along the line 
connecting the pulsar to the centroid of the X-ray emission in the 
torus, which corresponds roughly to the jet line (see Sect. 13. 2\ . 

Performing phase-resolved imaging (see Sect. 13. 4K we found 
that the on-pulse emission (where the Crab pulsar accounts for 
an important fraction of the flux) depends on energy: the on- 
pulse shift with respect to the pulsar position decreases by ~ 
30% between 20 and 200 keV. This effect can be easily explained 
by an increase in the pulsed fraction (see Sect. 13.3b . 

More interestingly, we find that during the off-pulse phase 
(where the emission from the pulsar is negligible) the centroid 
of the source is significantly offset from the position predicted 
from Chandra data (see Sect.|4]i. This indicates that the morphol- 
ogy of the source changes around the break energy. Consistently, 
we show that a cut-off in the electron spectrum in the outer re- 
gions of the X-ray nebula, caused by synchrotron cooling and i n 
agreement with the predictions of lAtovan & Aharonianl (fl996). 
can reproduce the observed shift. 

The centroid of the nebular emission around the break en- 
ergy coincides with the inner ring, which is int erpreted as the 
locati on to be the wind termination shock ( Weisskop f et al.l 
2000). Therefore, our data imply that above the break energy 
only the strong shock regions are responsible for the emission. 
Comparing our results with the predictions of a simple model for 
the spectral evolution of the nebula, we find that the electrons in 



the outer regions of the torus (d > 0.35 pc away from the pulsar) 
probably reach a maximal energy close to 10 14 eV. This result 
agree s with theoretical studies carried out using MHD simula- 
tions (IVobi et alJl2008l) . In the near future, because of a higher 
angular resolution by a factor ~ 20 NuSTAR will be able to re- 
solve the Crab pulsar/nebula complex up to ~ 80 keV, which will 
allow us to probe the electron spectrum in the different regions 
of the nebula close to the break energy. 
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